Subphotospheric Neutrinos from Gamma-Ray Bursts: The Role of Neutrons 
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Relativistic outflows with neutrons lead to hadronuclear reactions, and resulting subphotospheric 
7 rays may explain prompt emission of 7-ray bursts. In this scenario, quasi-thermal neutrinos in 
the 10-100 GeV range should be generated, and they may even have a high-energy tail formed 
by neutron-proton conversion or shock acceleration mechanisms. We demonstrate the importance 
of dedicated searches with DeepCore+IceCube, though such analyses have not been performed. 
Detections allow us to see roles of relativistic neutrons as well as effects of cosmic-ray acceleration. 
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I. INTRODUCTION 

Neutrinos provide an important messenger of explosive 
phenomena such as gamma-ray bursts (GRBs) and su- 
pernovae (SNc). They can probe physical processes deep 
under the photosphere (at which the Thomson opacity 
tt is unity), which cannot be seen by electromagnetic 
signals. Thermal MeV neutrinos are mainly radiated 
from proto-neutron stars or accretion disks at the core- 
collapse of massive stars, as detected from SN 1987A. On 
the other hand, assuming that outflows accompany ion 
acceleration via shocks or magnetic reconnections, high- 
energy neutrinos have theoretically been predicted above 
0.1 TeV drfll- They are detectable especially by neutrino 
detectors such as IceCube Q and KM3Net @. 

Many predictions on GRB neutrinos have been made 
for prompt 7-ray emission observed in the MeV range 
and its spectra can often be fitted by a smoothed bro- 
ken power law (the so-called Band function) [loj]. These 
neutrinos are also motivated by the hypothesis that 
GRBs arc sources of ultrahigh-energy cosmic rays (UHE- 
CRs) Q. Most of the works arc based on the classi- 
cal scenario, where 7 rays are attributed to optically- 
thin synchrotron emission from electrons accelerated at 
internal shocks [ll|. However, IceCube observations 
have constrained parameter space expected in this sce- 
nario [HI, [l]| , favoring large dissipation radii [l4| . More- 
over, the classical scenario has several difficulties, and the 
emission mechanism has been a long-standing issue [l5j |. 

One of the alternative scenarios is the photospheric 
scenario, where prom pt 7 ra y s are generated at subpho- 
tosphcrcs (tt > 1) [12l [l8[. This scenario has some 
advantages [l5j], an d observations have also indicated a 
thermal-like component [l9| . In particular, subphoto- 
spheric dissipation may originate from inelastic nucleon- 
neutron collisions just beyond the decoupling radius fl8j . 
This JVn collision model naturally leads to quasi-thermal 
7-ray spectra as a result of cascades and Coulomb heat- 
ing. Also, the neutron-loading in GRB jets is naturally 
expected in models of the GRB central engine [20j |. 

Importantly, quasi-thermal neutrinos are inevitably 



accompanied by those neutron-induced 7 rays, unlike 
conventional non-thermal neutrinos that rely on uncer- 
tain cosmic-ray acceleration mechanisms and may not be 
produced by radiation- mediated shocks 0, [HJ . They are 
expected at ~ 10 — 100 GeV, so not only IceCube [8| but 
also its low-energy extension DeepCore [21j are crucial for 
detections. In this work, we demonstrate that searching 
for the quasi-thermal neutrinos is relevant and detecting 
them allows us not only to test a neutron component but 
also to see if cosmic-ray acceleration occurs. We also dis- 
cuss the importance of low-luminosity or failed GRBs. 
We use Q x = Q/10 x in CGS units. 



II. QUASI-THERMAL EMISSIONS 

We consider neutron-loaded relativistic jets, which are 
naturally expected in GRB engines including accretion 
disks and proto-neutron stars [20(. In the baryonic fire- 
ball scenario, their final Lorentz factor achieves T r « rj, 
which is the initial value of random internal energy per 
particle. Initially, protons and neutrons are well coupled, 
and they should be decoupled at subphotospheres. If the 
decoupling happens before coasting, the neutron compo- 
nent makes the slower flow with T s = T n . Then, the 
faster flow overtakes the slower flow, leading to internal 
collisions. Even if the coasting is earlier, inhomogeneity 
in the jet leads to internal shocks at r ss 2T 2 s ri [ll[ , where 
Ti is the jet basis. 

Considering a collision between outflows with F r and 
r s , the Lorentz factor of the interaction shell is T sa 
^jT r T s / (r s /r r + Tp n ) when the faster flow collides with 



the 



Here 



neutron flow at r pn < 1 [181 . |22ij . i pn 

,„(r/r„) is the opacity for the pn reaction, where 



) and L n is the neutron luminos- 



ity. The pn cross section is a pn ~ 3 x 10 



" 26 cm 2 



which 



is comparable to the pp cross section 
completely merge, we have T 



When the flows 



vTVTV 

As relativistic neutrons loaded in the slower flow go 
through, the kinetic energy of the faster flow should 
be dissipated via inelastic nucleon-neutron collisions. 
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Then, mesons and muons arc produced, which d ecay 
into 7 rays, electrons (positrons) and neutrinos [23| . 
High-energy 7 rays cannot avoid the 77 process, and 
they induce cascades in the source. The cascades in- 
crease the number of pairs, which are heated by pro- 
tons via Coulomb collisions [18[ and/or by plasma/MHD 
waves [24| . Such so-called slow- heating models can also 
explain the spectra of the prompt emission. Since pairs 
are efficiently injected by the cascades, the Thomson 
opacity can be enhanced compared to for baryon- 
associated electrons. When the pair density is de- 
termined by the balance between Coulomb heating 
and inverse-Compton (IC) cooling, one obtains tt ~ 
23(y±/0.2) 1/2 (L/5L„) 1/2 (r/5r„)- 1 r p „ 0, where Y± is 
the pair yield and L is the kinetic luminosity of the flow 
with r. When 7 rays are produced at subphotospheres, 
the resulting emission should be modified by Compton 
scatterings before leaving the emission region, and its 
spectra can be compatible with observations. 

Protons in the faster flow are thcrmalizcd by pn colli- 
sions, leading to quasi-thermal nuclcons with rclativistic 
temperatures. Their typical energy is 
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K p r re i rripC 2 ~ 1.4 GeV T 
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(1) 



in the comoving frame. Here n p w 0.5 is the proton 
inelasticity and r rc i w 0.5(r/r s + r s /r). Collisions at 
1 generate efficient 7-ray emissions [l8j], so the 



' pn 

energy of quasi-thermal nucleons is ~ 0.5£]^°, where 
£jv° is the kinetic energy that dissipates. The significant 
fraction of the dissipated and trapped energy (~ 0.5£]^°) 
can be released as 7 rays. Assuming that half of the 
energy is used for adiabatic expansion, we may expect 
£^°/£^° ~ 0.25, though it can be smaller. 

Neutrinos easily leave the source, providing the more 
direct signal. When the faster flow is decelerated by col- 
lisions with neutrons, the observed neutrino energy is 
typically E v 0.05T r m p c 2 . The emission comes from 
the flow with V, so it is useful to use T and T' lel « 
0.5(r r /r + r/r r ) « r re iT p „. Using r r ^> r s , we have 



Ef « o.irr' 



rcl" t P<- 



150 GeV r 2 . 7 r; 
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(2) 



implying - 30 - 300 GeV neutrinos for T - 10 2 - 10 3 . A 
neutrino typically carries 1 /4 of the pion energy but this 
energy fraction ranges from to 0.43 in decay, and 
the high-energy tail is important for the detectability of 
neutrinos. It is reasonable to take r re ] ~ a few. 

Particle acceleration. — Quasi-thermal neutrino emis- 
sion does not rely on uncertain cosmic-ray acceleration 
mechanisms. One may also invoke shock acceleration 
when both flows contain protons. Note that neutrons 
can go through the faster flow when the neutron penetra- 
tion depth ~ (K p CT„ p njv) 1 is longer than r/T, so inter- 
nal shocks beyond the decoupling radius are considered. 
Cosmic-ray acceleration would be suppressed at shocks 
with tt 1 @, EH- However, collisionless (radiation- 
unmediatcd) shocks and ion acceleration become possi- 
ble especially at the reverse shock once tt < a few by, 



e.g., reducing the opacity enhancement with smaller Y±. 
We also consider this situation, assuming the proton ac- 
celeration time t acc w 1-KEpl (eBc) 0. The spectrum is 
expected to be a power- law (PL) with the index s, and 
its normalization is given by the efficiency e acc . 

Another plausible possibility is the neutron-proton 
conversion (NPC) mechanism [25|, [26|. Hadronuclcar re- 
actions with incoming neutrons inevitably generate rela- 
tivistic nuleons in the downstream. Protons are quickly 
isotropized by magnetic fields while experiencing the pp 
reaction. Then, as in the shock acceleration mechanism, 
a fraction of neutrons produced as a result of n — > p — > n 
can go back into the upstream, and they are overtaken 
by the shock front after the next conversion. The shock 



acceleration theory predicts an initial boost by 2r^ el [27 1 . 
Similarly, quasi-thermal protons can be boosted despite 
energy losses, and we have 



e^ c » o.05r(2<r4 1 ) 
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180 GeV r 2 .rr 3 cli0 . 5 . 



(3) 



In principle, further boosts could be relevant when r rc i is 
large. But using Eq. (3) is enough due to supprresion by 
other cooling processes. Note that the NPC acceleration 
time must be comparable to the hadronuclcar collision 
time. Obviously, the NPC mechanism is efficient only 
when Tp n is not small. Naive considerations lead to the 
efficiency e npc sa f/ npc (r^ el /8)min[l, rjf ] and Monte-Carlo 
simulations suggest <7 npa ~ 0.01 — 0.1 j26j |. 



III. NEUTRINO EMISSIONS 

We numerically obtain the neutrino spectra with 
Geant4 [2^], following For the quasi-thermal com- 
ponent, we calculate hadronuclear reactions between 
cold nucleons with T rc i and incoming neutrons with 
r re i. Quasi-thermal nucleons with are assumed to 
be isotropic and to lose their energies via various cooling 
processes. Here we consider energy losses by Coulomb 
collisions, hadronuclear reactions, Bethe-Hcitlcr process, 
photomeson production, synchrotron/IC emission and 
adiabatic expansion. Decaying pions, kaons and muons 
also lose their energies via hadronic processes, radiative 
cooling and adiabatic expansion, which are treated by 
solving kinetic equations. There are four principal pa- 
rameters, which are set to V — 600, r rc i = 3, r pn = 1 
and £%° = 45^° (= 2£^°). The former two determine 
neutrino energies, while the latter two do the fluence nor- 
malization. We also use L„ = 2 x 10 51 erg s _1 , where 



Tpn = 1 corresponds to r ~ 1.1 x 10 cm. The photon 

temperature is given by T m 530 eV L\ ^, r i\ ^2 7^ 2 anc ^ 
the magnetic field fraction is set to eb — 0.01 though they 
are not critical for results. 

We also consider the NPC component that is approx- 
imated by a mono-energetic distribution with 0.5r 2 el e^. 
This is enough for our purpose, and the normalization 
is set by e npc = 0.1. For the PL component, we de- 
termine the proton maximum energy £p lax by i acc < 
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FIG. 1: The energy fluence of + from a high- luminosity 
GRB (with 4 S ° = 10 53 - 5 erg) at z = 0.1. The ANB in 30 s is 
shown by the dot-dashed curve. 

min[i p , tdyn], where t p is the proton cooling time and 
idyn ~ r/(Tc) is the dynamical time @. Our parameters 
imply £™ ax ~ 10 5 ' 6 GeV. Motivated by recent numerical 
simulations (2^, we adopt s = 2.1 and e acc = 0.3. 

The results for a high-luminosity GRB at z = 0.1 are 
shown in Fig. 1. As expected in Eq. (1), quasi-thermal 
neutrinos have a peak at ~ 100 GeV. The NPC compo- 
nent enhances a high-energy tail, but it is not very rele- 
vant for our conservative value of r re i. The PL compo- 
nent is prominent above TeV, and hadronuclear reactions 
give a dominant contribution especially for steeper spec- 
tral indices. The photomeson production is also quite 
efficient, but the fluence is largely suppressed by strong 
cooling of mesons and muons. We also show the atmo- 
spheric neutrino background (ANB) (3(| assuming that 
the angular window of max[0 2 , 7r# 2 ], with = 1 deg and 
the kinematic angle Q v « 1.5 deg ^/TcV jE v . 

Detecting neutrinos from one GRB requires nearby 
bursts. But most of these are much less energetic bursts 
like GRB 060218 which may originate from low T jets 
or shock breakout from jet-driven SNe [3l|. Note that 
hadronuclear collisions may occur even inside the stellar 
envelope, so subphotospheric neutrinos are expected from 
choked jets [H, HH as well as successful jets. The results 
for a low-luminosity GRB at D = 10 Mpc are shown 
in Fig. 2, with T = 30, r ro i = 5, and a sub-parameter 
L n = 2 x 10 46 erg s _1 . Quasi-thermal neutrinos are ex- 
pected around 10 GeV, which also demonstrates lower T 
cases. The NPC component, which is prominent above 
100 GeV due to higher r re i, is shown with e npc = 0.3. 



IV. NEUTRINO DETECTION 

Since IceCube is not sensitive at 10 — 100 GeV, 
including DeepCore is essential to see quasi-thermal 
neutrinos. The neutrino effective area of Deep- 
Corc+IccCube at 10 - 100 GeV is roughly w 
10 15 cm 2 (£„/100 GeV) 2 [U, so detections at E v rc- 



FIG. 2: The same as Fig. 1, but for a low-luminosity GRB 
(with 4 S ° = 10 50 erg) at D = 10 Mpc. The ANB in 1000 s is 
shown by the dot-dashed curve. 

quire El4> v > 5 x 10" 3 erg enr 2 {E U /10Q GeV) -1 . Only 
energetic and nearby GRBs can be seen, and a few events 
are detectable in the case shown in Fig. 1. 

Hence, it is critical to make dedicated stacking anal- 
yses for GRBs detected by 7-ray satellites. Although 
such analyses have been done around PeV energies for 
the classical scenario [lU, [l3| , but not at < 1 TeV for the 
photosphcric scenario. To demonstrate how to search for 
subphotospheric neutrinos, we use the fluence distribu- 
tion obtained by Fermi-GBM (see Fig. 7 in |34j|). GBM 
detected 400 long bursts in two years, and we assume 
that GBM sees 2000 bursts in the northern hemisphere 
in 20 years. To discover the signal, the number of events 
has to be enough and the signal-to-background should be 
sufficiently large. From Fig. 1, the ANB at - 100 GeV 
is ~ 10 -6 erg cm -2 , so the fluence threshold for stack- 
ing should be > 10~ 6 erg cm~ 2 . Taking thresholds of 
< 10~ 6 erg cm -2 is not useful since the integral fluence 
distribution is flat there, while using higher thresholds 
is not very essential since the smaller number of more 
energetic bursts is compensated by higher fluences. 

How we normalize the fluence is crucial. In the classical 
scenario, the normalization is given by the GRB-UHECR 
hypothesis [l[ or a cosmic-ray loading parameter 0]. In 
this work, analogously to the hadronic model for an extra 
GeV component (35[, we use the observed 7-ray fluence 
as E^4> 1 (x since subphotospheric 7 rays are assumed 
to be responsible for the prompt emission. Second, the 
meson production efficiency / P7 affects the fluence. In 
the classical scenario, / P7 is sensitive to r and T that 
are uncertain I3J. In our model, dissipation should oc- 
cur at tt ~ 1 [l||, and efficient 7-ray production should 
accompany neutrinos. Finally, the typical neutrino en- 
ergy depends on uncertain T and z. For simplicity, we 
fix r = 600 and z = 1. Similar assumptions were also 
made in analyses for the classical scenario [l2|, [l3| , where 
the typical energy depends on T and r (for sufficiently 
high f pl ) as well as L 7 and break energy @,||- 

The expected number N of detected v ^ + events is 
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FIG. 3: The number of detected events for v^+v^. Coincident 
20 yr observations with DeepCore+IceCube and GBM are 
assumed. The dot-dashed curve is the ANB. 
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FIG. 4: The same as Fig. 3 but for v R + v R that are observed 
via cascades. The angular resolution for the ANB is assumed 
to be 20 deg. The dot-dashed curve is the ANB. 

shown in Fig. 3, with the threshold 10 -5 5 erg cm~ 2 . The 
effective areas of DeepCore and IceCube are taken from 
pl| and @, respectively. We predict that a few events 
will be detected by analyzing ~ 1000 GRBs stacked in 
20 yr observations. Given the fluence at E v , N roughly 
decreases with E v in DeepCore, so stacking ~ 5000 GRBs 
is needed to find quasi-thermal neutrinos if all GRBs have 
r ~ 100. On the other hand, we can expect higher T for 
energetic bursts, as suggested in LAT GRBs (3g|. Quasi- 
thermal neutrinos lead to plateaus below ~ 100 GeV due 
to their narrow distribution, and they are more promising 
than non-thermal neutrinos for subphotospheric emis- 
sion. Stacking ~ 2000 bursts allows us to see a PL com- 
ponent without depending on T. 

Muon neutrinos are mainly detected from muon tracks, 
whereas electron neutrinos are seen via cascades (Fig. 4). 



The ANB is more severe since the angular resolution is 
worse. But better reconstruction techniques can improve 
the detectability significantly, e.g., if the low-energy ex- 
tension of KM3Net could achieve ~ 5 deg [37j • 

V. DISCUSSIONS AND IMPLICATIONS 

If neutrons play a major role in generating prompt 
7 rays, these should be accompanied by quasi-thermal 
10 — 100 GeV neutrinos. The signal is much more ro- 
bust than more conventional non-thermal neutrinos that 
rely on uncertain cosmic-ray acceleration mechanisms. 
Dedicated searches have not been done and using only 
IceCube is insufficient. We strongly encourage stack- 
ing analyses with low-energy extensions of IceCube and 
KM3Net, and detections are possible in decades with 
DeepCore- like detectors. Nearby low-luminosity and/or 
failed GRBs can also be interesting targets for revealing 
the GRB-SN connection. 

Neutrons play various roles [HI, including dissipation 
via pn collisions and production of quasi-thermal parti- 
cles. These rclativistic particles may naturally become 
seeds injected into the cosmic-ray acceleration processes. 
In addition, neutrons may generate magnetic fields via 
np conversions. As neutrons go through the unmagne- 
tized faster flow, they inject proton beams and quasi- 
thermal protons with relativistic temperatures. In par- 
ticular, plasma anisotropics may lead to filamentation or 
Weibel instabilities, making the faster flow magnetized. 
The magnetic fields are important for scattering of parti- 
cles [24[ as well as synchrotron emission of electrons [l8[ . 
Detections of ~ 10 — 100 GeV neutrinos would provide 
hints of these effects of neutrons. At least, they would 
provide evidence of subphotospheric dissipation leading 
to quasi-thermal nuclcons. Detecting > 0.1 TeV neutri- 
nos could provide insights into the NPC mechanism and 
cosmic-ray acceleration under the photosphere. 

Note that our results are consistent with IceCube non- 
detections of neutrinos even though the Nn collision 
model requires large kinetic luminosities at the colli- 
sion radii. If the jet is initially Poynting-dominated, for 
this model to work, the magnetic energy must be some- 
how quickly converted. If not, neutron-induced emissions 
should be less, and other types of magnetic dissipation 
would be required for the prompt emission J3j|. On the 
other hand, if ions are efficiently accelerated via magnetic 
dissipation, we may expect the non-thermal neutrinos. 
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